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Renormalization constants (Z-factors ) of vector and axial-vector currents are determined non- 
perturbatively in quenched QCD for a renormalization group improved gauge action and a tadpole 
improved clover quark action using the Schrodinger functional method. Non-perturbative values of 
Z-factors turn out to be smaller than one- loop perturbative values by 0(15%) at lattice spacing 
' of a^ 1 pa 1 GeV. The pseudoscalar and vector meson decay constants calculated with the non- 

perturbative Z-factors show a much better scaling behavior compared to previous results obtained 
£SJ ■ with tadpole improved one-loop Z-factors. In particular, the non-perturbative Z-factors normalized 

^ ' at infinite physical volume show that scaling violation of the decay constants are within about 10% 

up to the lattice spacing a" 1 ~ 1 GeV. The continuum estimates obtained from data in the range 
a" 1 = 1 - 2 GeV agree with those determined from finer lattices (a 1 ~ 2 — 4 GeV) with the 
standard action. 
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I. INTRODUCTION 



Reliable prediction of physical quantities from lattice QCD calculations requires a good control of scaling violations. 
For this purpose, several improved actions have been tested and applied to large scale systematic simulations. For 
most of physical quantities, such as quark masses and hadronic matrix elements, one, in addition, has to calculate 
renormalization constants (Z- factors) . Non-perturbative methods to determine various Z- factors have been developed 
and utilized for several actions. 

The CP-PACS collaboration carried out systematic simulations for two-flavor full QCD [1] using a renormalization 
group (RG) improved gauge action [2] and a tadpole improved [3] clover quark action [4] and reported results for 
the spectrum, quark masses and meson decay constants. Since non-perturbative Z-factors were not available for 
this action combination, analysis had to rely on one-loop perturbative values. The result showed that the meson 
decay constants in both quenched and full QCD suffer from large scaling violations at a -1 ^1 — 2 GeV which 
hinder continuum extrapolations. It was not clear if this observation hinted at an inherent difficulty of improved 
actions for matrix elements at coarse lattice spacings or a perturbative treatment of the Z-factors was the issue. A 
non-perturbative determination of Z-factors was evidently needed. 

As a first step toward a systematic study of non-perturbative renormalization for this action, we apply the 
Schrodinger functional (SF) method [5-8] to calculations of Z-factors for vector (Zy) and axial-vector (Za) cur- 
rents in quenched QCD with the same improved action. The SF method has been applied to the non-perturbatively 
0(a) improved Wilson action. On the other hand, our action combination has 0(ag ) error, since the coefficient of 
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the clover term is determined by a tadpole improved perturbation theory to one-loop order. Therefore our study 
involves an examination whether the SF method successfully works for this action 

In the SF method, Z-factors are determined at various couplings for a fixed physical size L. In this case, the 
Z-factors contain terms of 0(ag 4 /L) in addition to that of O(aA). We expect that one can remove the 0(ag 4 /L) 
terms by taking the infinite volume limit. We calculate Z- factors both at a finite fixed physical volume and at the 
infinite volume and compare scaling properties of decay constants with these two choices of Z-factors. 

When we calculate Z-factors for large couplings and on large lattices, we encounter "exceptional configurations" 
for which observables take abnormally large values. We estimate systematic uncertainties due to such configurations, 
propagate them to systematic error estimations of decay constants, and argue that they do not alter our conclusions 
on scaling properties. 

The organization of this paper is as follows. In Sec. II we describe calculational method focusing on features of 
applying the SF method to the RG improved action. Detail of analysis for Zy and Za are given in Sec. III. With 
these ^-factors, we study scaling behavior of decay constants of vector (/ p ) and pseudoscalar (f v ) mesons in Sec. IV. 
Sec. V is devoted to conclusions. Parts of this work have already been reported in Refs. [9] and [10]. 



II. CALCULATIONAL METHOD 



A. Action and currents 



The RG improved gauge action we employ has the form 



K x,fj.,i^ x,fj,,v ) 



S g = - { c 2^ Up^ v {x) + ci 2^ U R .^(x) } , ( 1 ) 

where (3 = 6/g 2 with the bare coupling constant g, and Up (Up) is the trace of the product of link variables around a 
plaquette (6-link rectangular loop). In a sum over loops, each oriented loop appears once. The coefficients c\ = —0.331 
and c = 1 — 8ci = 3.648 are fixed by an approximate RG analysis [2]. 
The clover quark action [4] is defined by 

S, = 2?(aOJW(l>). (2) 

x,y 

D X ,y = S x .y - K^{(1 - J^Ux.^Sx + ft.y + {^-+l l j)Ul^5 X ^ + lx} - 6 X ,y C S W « ^ ^ ^ ' ^ 

where n is the hopping parameter and F^ v is a lattice discretization of field strength. For the clover coefficient csw, 
we adopt a tadpole improved value with one-loop estimate of the plaquette (W lxl ) — -| (Up^ v ) given by 

csw = <VF 1X1 )" 3/4 = (1 - 0.8412//3)- 3 / 4 , (4) 

since the one-loop value of the plaquette reproduces measured values well. 

For Zy, we investigate vector Ward identity of the unimproved current V® since the particular SF setup we use 
yields the same value of Zy for improved and unimproved currents. On the other hand, we study an improved current 
for Za defined by 



(Aj)^x) = A%{x) + ac A -(d; + d,)P a (x), (5) 



where d* and are lattice forward and backward derivatives and P a (x) is the pseudoscalar density. We use one-loop 
value for the improvement coefficient ca — — 0.0038g 2 [11]. 



B. Implementation for RG improved action 

We follow the SF method of Ref [8] for the RG improved gauge action. Taking an L 3 x T lattice, we impose the 
periodic boundary condition in the spatial directions and the Dirichlet boundary condition in the temporal direction. 
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The boundary counter terms for the RG improved action are determined [12] so that the classical field equations 
are satisfied. The total action reads 



S s=qzZzZ \ co^w»Axi)Up,n,v{x) +c 1 ^2,w* l/ (x A )U Rttlv {x) \ , (6) 

x x 4 =0 I p*.y ii,v ) 

with the weight factors w^ v {x/C) and w pu (x4) given by 

c t (X4 = or X4 = T — a) and (fj, = 4 or v = 4) 

. x 4 = T and {pi = 4 or v = 4) 

'"' ' ^c s (a* = or X4, = T) and (/j + 4 and i/ ^ 4) ' {l ' 

1 otherwise. 

| when 2 links of the rectangular touch a boundary 
wjf„(x4,) = i when the rectangular is completely included in a boundary (8) 
1 otherwise. 

We take the tree- level value of c s = c t = 1. (One-loop values of c s and c t [13] were not known when we started this 
work.) Boundary counter terms are not included for the clover quark action. In other words, the boundary coefficients 
c s and Ct of Rcf. [8] are set to the tree-level values, c s — c t = 1. 



C. Details of calculational method 



We calculate Z- factors for the range (3 = 8.0 — 2.2. The smallest value (3 — 2.2 is chosen so as to cover the range 
(3 = 2.575 — 2.184 (inverse lattice spacing a -1 w 1 — 2 GeV) where data for f p and f„ exist [1]. Simulations are made 
with a 5-hit pseudo- heat-bath algorithm mixed with an over-relaxation algorithm in the ratio of 1:4. We analyze 200 
- 20000 configurations separated by 100 sweeps each. Lattice geometry is set to T — 2L for both Zy and Za- At 
each (3, Z-factors are determined for at least two lattice sizes in order to interpolate or extrapolate them to a fixed 
physical volume. Lattice size and number of analyzed configurations are listed in Table I. 

Calculations of Z-factors are carried out at zero quark mass, where the quark mass m q is defined by the PCAC 
relation d^A^ = 2m q P a in the continuum notation. The actual procedure to measure m q is the same as in Ref. [7]. 
We define a time-dependent quark mass by 



| (<9o + d ) fA(xi) + c A ad^d f P (x 4 ) 



m q ( Xi ) = ^ " of (r] > ( fl ) 

where Ja is the axial-vector and pseudoscalar correlator and fp is the pseudoscalar and pseudoscalar correlator. 
Quark mass m q is defined by an average of 77^(24) over a range of time slice around X4 — T/2; 



1 an ( T \ 

— zZ m " 2 +t )' (10) 



H 2n + 

where n defines the range and depends on simulation parameters. 

Z- factors are determined from correlators of pseudoscalar operators P° and P' a at the boundaries and/or currents 
V® and (Ai)® at a finite time slice. We employ the notations defined by 

h = -^(P' a P a ), (11) 

fv{Xi) ^E^^w^), (12) 

X 

fAA(x4, V4) = eabC " Cde ( P ' d ( A l) a (x)(Al) b (y)P e ) • (13) 

Renormalization constants are then extracted from 
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2v = TTr-y ( 14 ) 

Za = \ ~f ( 3T 5T\ ' ( 15 ) 

V Jaa{— ,— ) 

Lattice spacing necessary to set the physical size is determined from the string tension = 440 MeV. We fit 
values of a\fa [1,14] to a fitting form [15] 

(aV^)(/3) = /(/3)(l + C2 a 2 (/3)+ C4 a 4 (/3) + ---)/co, a(/J) = (16) 

/(.Pi) 

where /(/?) is the two-loop scaling function of the SU(3) gauge theory, c„'s are parameters to describe deviation from 
the two-loop scaling, and p\ is a reference point. Choosing (3\ = 2.4, the parameters 

c = 0.5443(97), c 2 = 0.390(38), c 4 = 0.049(12) (17) 

reproduce measured values well. 



III. RENORMALIZATION CONSTANTS 



A. Results at simulation points 

Calculations of Z-factors do not present any difficulty from (3 = 8.0 down to 2.6 for all lattice sizes, and at f3 = 2.4 
and 2.2 on a small lattice of 4 3 x 8. For a larger 8 3 x 16 lattice at /? = 2.4 and 2.2, however, we encounter the issue of 
"exceptional configurations" . Deferring discussions of this issue, let us first summarize results for the non-exceptional 
case. 

We first determine k c for each (3 and L/a. For this purpose simulations are carried out at several values of k around 
an estimated k c in which m q is determined by eq.(10); we employ the fitting range given by n = for L/a = 4 lattices, 
n=lfor£/a = 8 lattices at (3 = 2.4 and 2.2, and n = 2 for L/a = 8 at other values of (3 and L/a — 16 lattices. We 
determine k c by a linear fit in n as illustrated in Fig. 1 for an 8 3 x 16 lattice at (3 = 2.6. 

Results for n c are given in Table I. When we perform simulations for Z-factors at k c , m q is also calculated for 
confirmation. Typical result for the effective quark mass m q (x4) is shown in Fig. 2. It exhibits a reasonable plateau. 
Values of m q at n c are given in the fourth column of Table I. They are consistent with zero within an accuracy of 
10 -3 at worst. Since k c is tuned well, errors of k c are not taken into account in the error estimation of Z-factors. 

Figures 3 and 4 show typical plots of Z v and Z4 as a function of time slice. We observe good plateaux for all 
cases, except for a small temporal lattice size of L/a — 8 at strong coupling of (3 = 2.4 and 2.2. We thus find that 
the SF method works successfully to our action, albeit there are 0(ag 4 ) errors rather than 0(a 2 ) errors. Results for 
Z-factors at simulation points are given in Table I. 

We now discuss the issue of "exceptional configurations". Among our parameter sets, anomalously large values 
appear in the ensemble of hadron correlators on an 8 3 x 16 lattice at (3 = 2.4 and 2.2. We illustrate the situation in 
Fig. 5, where we plot the time history of j 'a and /p at /3 = 2.4. No large spikes appear for larger (3 on the same size 
lattice, as shown in Fig. 6 for data at (3 = 2.6. Note that a much finer vertical scale is employed in the latter figure. 
Large spikes appear also in the ensembles of /1, fy and j ' aa necessary for evaluating Z- factors. Such Exceptional 
configurations make it difficult to determine m qi k c and Z- factors precisely. 

We suspect that the spikes are caused by the appearance of very small or even negative eigenvalues of the Dirac 
operator toward strong coupling. Such eigenvalues would be suppressed in full QCD by the quark determinant, and 
in this sense we expect exceptional configurations to be an artifact of quenched QCD. In quenched QCD, however, 
"exceptional configurations" cannot be distinguished from "normal" ones on some rigorous basis. In fact, histograms 
of /'s have a long tail toward very large values as shown in Fig. 7. We then adopt the strategy of estimating Z-factors 
by removing from the ensemble average configurations having values above some cutoff. The uncertainties associated 
with this procedure are estimated by varying the cutoff, and will be propagated to systematic errors of Z-factors at 
fixed physical volume. Detailed description of this procedure will be given in Sec. IIIC. 
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B. Z- factors at fixed physical volume 



We plot Zy and Za determined for various sizes and in Figs. 8 and 9, respectively, as a function of a/L. There 
are three or four data points at = 8.0 and 2.8 for Zy, and at = 2.8 for Za- These data show a linear behavior 
in a/L, which is consistent with the expectation that Z-factors for our action have 0{ag / L) errors. Therefore we 
adopt a linear ansatz to extrapolate or interpolate Z-factors to the physical volume of L — 0.8 fm (normalized at 
= 2.6 and L/a = 8) and to L = oo, as shown in these figures. We denote Z-factors at a fixed physical volume as 
Zy F A L ~°' Stm and Zy^ i_0 °. Numerical values are given in Table II. 

In Fig. 10 we plot Zy as a function of g 2 . Making a Pade fit, we obtain 

7 SF,L=o. 8 fm _ l- 0.302225g 2 + 0.011034g 4 

Zy ~ 1 - 0.23943V ' { *> 

7 sf.l=oo _ 1 - 0.365802g 2 + 0.015016. 9 4 
v 1 - 0.303008. 9 2 ' [ ' 

where we have imposed a constraint that the Pade fits reproduce the one-loop perturbative result Zy = 1 — 
0.062794 3 2 [16] up to 0(g 2 ). 

The range of coupling where there are data for the vector meson decay constant is marked by the two vertical 
dashed lines in Fig. 10. Over this range, Z V F ' L ~°° becomes increasingly smaller (by about 8 - 21%) compared to the 
one-loop value (dashed line) and the tadpole- improved value (crosses). We also observe that Zy exhibits a sizable 
volume dependence toward strong coupling, e.g., for < 2.8. This will have an important consequence on the scaling 
property of f p as discussed in Sec. IV. 

Another method to estimate Zy non-perturbatively [17] utilizes the non-renormalizability of the conserved vector 
current 

Vf(x) = ^{^{x + a(i)Ul^{ li + l)i ) {x)+^{x)U x ^{ li -l)ij{x + aix)}, (20) 



and defines 



7 NPC_ „„ Ex(0|^ C (x,Z4)^|0> 



V £x«W(x,Z4W|0> • (21) 

Results for Zy FC obtained for our action combination [1] are also overlaid in Fig. 10 (open diamonds). They are 
much smaller than those from the SF method. We interpret that the large difference originates from large 0(a) errors 
in Zy P . From the viewpoint of 0{a) improvement of operators, divergence of a tensor operator d^T^ should be 
added to both the local current V$ and the conserved current [18]. However, the improvement operator djToj, 
necessary for Vq in the SF scheme, automatically drops out since it is a spatial total divergence. In other words, Zy 
from the SF method is 0(a) improved, whereas Zy PC is not. 
In Fig. 11, we plot results for Za and Pade fits which read 



~ l- 0.220946 3 2 ' 1 ' 

7 sf.l =0 o _ l- 0.334232g 2 + 0.011710/ 



In the Pade fits, we use a constraint from one-loop perturbation theory that Za — 1 — 0.056630.g 2 + 0{g i ) [16]. We 
observe that Z^' L ~°° for the range we have data for / w are smaller than the tadpole improved value by 6 - 14%. 



C. Systematic error from exceptional configurations 

Exceptional configurations affect Z- factors in two way, firstly by changing the value of m q and hence that of n c , 
and secondly by directly affecting the value of Z-factors themselves. Hence, in order to estimate the uncertainties 
of Z-factors due to exceptional configurations at large couplings, we investigate how m q and Z-factors change if we 
discard configurations having values of relevant correlators larger than some cutoff. 

Fi gure 12 shows this test for m q at — 2.4 for which a cutoff is set for the value of fp. As the cutoff fp^t 
is increased, m q gradually decreases, becomes almost stable around /p.cut = 300 and then the error of m q becomes 
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large. With this observation in mind, we have estimated k c as the point where m q for /pcut = 300 is consistent 
with zero, as shown in Fig. 12. The uncertainty in m q at n c is estimated by varying f Pcu t fr° m 200 to 1000 and 
turns out to be +0.00028 > m q > —0.00182. The same procedure at (3 = 2.2 with / pcu ^=500 gives the uncertainty 
+0.00344 > m q > -0.00170. We note that the number of configurations discarded is 37 (191) of the total of 10000 
(20000) configurations at /? = 2.4 (2.2). 

The uncertainty of m q is translated to uncertainties of Z-factors. To do this, we carry out two additional simulations 
at k's slightly above and below n c . Figures 13 and 14 show how Z-factors depend on m q ; Zy is very stable 
against variation of m q , while Za shows a more prominent dependence. We fit the m q dependence of Z- factors 
by a linear function, and uncertainties of Z-factors are estimated by the difference of the central value and the 
maximum/minimum value for the range of error of m q . The uncertainties are given in Table III under the column 
8Z m " . 

We also estimate uncertainties in the statistical averaging of Z-factors themselves by varying the cutoff of f\. 
Figure 15 shows the ensemble average of f\ and fy together with Zy as a function of the cutoff .fi. cut . Though both 
fi and fy increase as /i, cu t, their ratio Zy is very stable around /i iCUt = 5.0 reflecting the fact that /i and fy arc 
highly and positively correlated. We determine the central value from /i. cu t = 5.0 and estimate errors by varying 
A, cut = 2.0 to 10.0. Figure 16 shows a similar test for Za- The cutoff dependence is more conspicuous than for 
Zy. Uncertainties thus estimated are listed in Table III under the column SZ except \ Note that we discard 47 (221) 
configurations at (3 = 2.4 (2.2). 

The two uncertainties SZ m " and $Z except - are simply added to estimate the total uncertainty SZ for L/a = 8 lattices. 
We then propagate them to uncertainties of Z-factors at fixed physical volume, listed in Table III. Uncertainties at 
L = 0.8 fm are smaller than those on L/a = 8 lattices, because the physical size is located between L/a = 8 and 
L/a = 4 lattices. Uncertainties of Z- factors at L = oo, enlarged by extrapolations, are larger than the statistical error 

£ ystat^L— oo 



D. Z- factors at simulation points for meson decay constants 

For discussions of scaling properties of meson decay constants, we need Z-factors at (3 values where raw data of 
/„ and A are taken. We evaluate them using the Pade fits obtained in sec. IIIB together with the estimates of 
uncertainties from exceptional configurations. 

For the latter purpose, we repeat Pade fits varying the Z-factor at j3 = 2.4 within the range of its uncertainty, 
calculate the Z-factor at a target (3 value, and take the difference between this value and the central value as an 
estimate of the systematic error from exceptional configurations at (3 = 2.4. The systematic error from those at 
f3 = 2.2 is estimated similarly and added linearly to obtain the total systematic error. 

Our final results for Z-factors at the simulation points for meson decay constants are listed in Table. IV. The 
systematic errors in Z v F ' L ~ co and Z S ^' L ~°° are comparable to the statistical ones and within at most 2cr. The ratio 
of systematic errors to statistical ones in the table is slightly smaller than the ratio at simulation points for Z-factors, 
(3 = 2.4 and 2.2, because the Pade fits are stabilized with data at (3 > 2.6 where no exceptional configurations appear. 
In the table, we reproduce the one-loop tadpole improved values Zy P A together with Zy PC for comparison. 



IV. SCALING PROPERTIES OF MESON DECAY CONSTANTS 



We investigate the scaling property of / p , combining the unrenormalized values of f p [1] and four different choices 
of Zy, 1) one-loop tadpole improved value, Zy P , 2) values determined by the SF method for L = 0.8 fm, Z V F ' L ~ ' 8 m , 
3) those at L = oo, z£ F ' i= °°, and 4) Z£J PC in eq. (21) defined by the conserved current. Unrenormalized f p was 
determined in Ref. [1] for unimproved current. Hence, renormalized f p has 0(ag 2 ) error. 

Numerical values of the renormalized f p are listed in Table V. Statistical errors in the unrenormalized f p and in 
Zy F,L_0 ' 8 m or Zy F ' L ~°° are added in quadrature, while systematic uncertainties for the latter are simply propagated 
to that of renormalized f p . Fortunately, systematic uncertainties from exceptional configurations are sufficiently small. 
They are at most 0.4tr of statistical errors and do not affect the conclusions below. Therefore we ignore systematic 
uncertainties in the following discussions. 

Figure 17 presents the results for f p as a function of a. Open circles are the previous result obtained with tadpole 
improved one-loop Z-factor, which exhibits a sizable scaling violation. If instead one uses the non-perturbative Zy 
for L = oo, we observe a much better scaling behavior (filled triangles). The Z-factor evaluated for L = 0.8 fm lies in 
between the two results (stars). 
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We find this result to be very encouraging; it shows that an apparent large scaling violation seen with the use of 
Zy P is largely due to neglect of higher order contributions in the Z-factor, and that inherent 0(a) errors in the vector 
decay constant are small even up to a large lattice spacing of a -1 <~ 1 GeV. 

We extrapolate f p with Zy F ' L ~°° to the continuum limit linearly in a. The value in the continuum limit f p = 
201.7(2.0) MeV turns out to be consistent with that f p = 205.7(6.6) MeV from a high precision simulation with the 
standard action [19]. 

We note that a relatively large scaling violation is observed for f p with Zy PC (open diamonds). Since the f p with 
Zy PC is nothing but decay constant determined from the conserved current in eq. (20), we suspect that 0(a) large 
scaling violation exists in the conserved current. 

The scaling property of f„ is investigated in a similar manner. We note that unrenormalized f„ in Ref. [1] was 
determined for the improved current using a tadpole improved one-loop value of ca with MS coupling, whereas Z- 
factors in this work employ ca with bare coupling. Difference in renormalized f„ arising from the difference of ca 
remains 0(ag A ). 

Figure 18 shows f v versus a. In this figure, we have not included the systematic uncertainty since, being smaller 
than statistical one, it does not alter the conclusions below. Numerical values are listed in Table V. The conclusions 
on scaling behavior are similar to those of f p . Employing which includes terms to all orders in g 2 and does 

not have 0(ag 4 /L) error, scaling behavior turns out to be better than previous results with tadpole improved one-loop 
Z-factor indicated. Making a linear continuum extrapolation, we obtain a value / ff = 123.3(8.7) MeV consistent with 
that /„■ — 120.0(5.7) MeV calculated with the standard action at significantly weaker couplings (open squares). 

V. CONCLUSIONS 

The Schrodinger functional method has been applied to calculations of vector and axial-vector renormalization 
constants for the combination of a RG-improved gauge action and a tadpole-improved clover quark action. With the 
Z-factors determined non-perturbatively, an apparent large scaling violation in the range of lattice spacing a -1 = 
1 — 2 GeV in the meson decay constants previously observed with one-loop perturbative Z-factors is significantly 
reduced. We conclude that the improvement attempted with the gluon and quark actions employed in the present 
work is effective for the meson decay constants as well. 

We find that scaling of decay constants is best when one uses Z-factors normalized at infinite volume. This suggests 
that removing 0(ag 4 /L) error in Z-factors by the limiting procedure L — > oo is important for actions with 0(ag 4 ) 
error. 

The non-perturbative Z-factors have enabled us to determine values in the continuum limit of decay constants. We 
may expect that other hadronic matrix elements are also reliably extracted from lattice spacings much coarser than 
a -1 w 2 GeV for our action combination if one uses Z- factors determined by the Schrodinger functional method. 

Finally we recall that a large scaling violation for the meson decay constants has been observed also in two-flavor full 
QCD with the same action combination of quenched QCD considered in this article. Non-perturbative determination 
of Z-factors will therefore be interesting to pursue for this case. Exceptional configurations are expected to be absent 
in full QCD. Therefore the Schrodinger functional calculation would be more straightforward. Work in this direction 
is in progress. 
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TABLE I. Simulation parameters and Z-factors at simulation points. For 8 3 x 16 lattices at (3 = 2.2 and 2.4 where we 
encounter exceptional configurations, we quote only central values and statistical errors. No calculations of Za were made for 
16 3 x 32 lattice at (3 = 2.8 and 12 3 x 24 lattice at f3 = 8.0. 
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L 3 x T 




m q 


Z v 


#conf for Zv 


Za 


#conf for Za 


2.2 


4 3 x 8 


0.139281(56) 


-0.000529(495) 


0.7499(11) 


2000 


0.7667(52) 


2000 




8 3 x 16 


0.140570(15) 


0.001493(399) 


0.7067(12) 


20000 


0.7350(610) 


20000 


2.4 


4 3 x 8 


0.136933(21) 


-0.000764(337) 


0.7995(10) 


1000 


0.8150(24) 


2000 




8 3 x 16 


0.137481(04) 


-0.000130(395) 


0.7652(04) 


10000 


0.7525(204) 


10000 


2.6 


4 3 x 8 


0.135558(08) 


-0.000276(341) 


0.8265(05) 


2000 


0.8474(16) 


2000 




8 3 x 16 


0.135701(10) 


-0.000464(154) 


0.8056(05) 


500 


0.8208(18) 


2000 


2.8 


4 3 x 8 


0.134532(10) 


0.000002(396) 


0.8482(06) 


1000 


0.8667(18) 


1000 




8 3 x 16 


0.134515(07) 


-0.000105(139) 


0.8312(04) 


500 


0.8482(13) 


1000 




12 3 x 24 


0.134554(08) 


0.000525(111) 


0.8268(03) 


500 


0.8520(17) 


1000 




16 3 x 32 


0.134587(09) 


0.000663(139) 


0.8231(06) 


500 






3.125 


12 3 x 24 


0.133209(01) 


-0.000039(091) 


0.8540(02) 


500 


0.8689(09) 


500 




16 3 x 32 


0.133219(05) 


-0.000092(060) 


0.8527(01) 


500 


0.8672(08) 


500 


4.0 


8 3 x 16 


0.131094(04) 


-0.000302(086) 


0.8972(02) 


500 


0.9068(06) 


500 




16 3 x 32 


0.131083(01) 


0.000160(044) 


0.8938(01) 


500 


0.9060(04) 


300 


6.0 


8 3 x 16 


0.128898(03) 


-0.000191(061) 


0.9366(01) 


300 


0.9425(03) 


300 




16 3 x 32 


0.128891(01) 


-0.000060(023) 


0.9337(01) 


300 


0.9398(02) 


300 


8.0 


8 3 x 16 


0.127869(02) 


0.000138(041) 


0.9545(01) 


300 


0.9600(02) 


300 




12 3 x 24 


0.127870(01) 


0.000576(331) 


0.9523(01) 


300 








16 3 x 32 


0.127870(01) 


0.000021(017) 


0.95151(3) 


300 


0.9565(01) 


200 



TABLE II. Z-factors at fixed physical volumes. 

■SF,L=0.8fm ySF,L=oo 7 5F,L=0.8fm ySF,L=oo 



2.2 0.7495(15) 0.6635(26) 0.7664(52) 0.7033(1221) 

2.4 0.7783(05) 0.7309(13) 0.7763(127) 0.6900(409) 

2.6 0.8056(05) 0.7847(11) 0.8208(18) 0.7942(39) 

2.8 0.8279(02) 0.8180(12) 0.8511(13) 0.8386(22) 

3.125 0.8526(01) 0.8488(07) 0.8671(08) 0.8621(42) 

4.0 0.8917(02) 0.8904(03) 0.9055(08) 0.9052(10) 

6.0 0.9309(02) 0.9308(02) 0.9372(05) 0.9371(05) 

8.0 0.9485(01) 0.9485(01) 0.9530(03) 0.9530(03) 



TABLE III. Breakdown of systematic uncertainties in ^-factors. See text for details. 







5Z m " 




c ^except. 




8Z 




gr?SF,L=Q.8tm 


SZ SF,L=. 






Zv(/3 = 


2.2) 


+0.04%, 


-0.02% 


+0.31%,- 


-0.06% 


+0.35%,- 


-0.08% 


+0.0013%, -0.0058% 


+0.75%,- 


-0.17% 


0.39% 


Z V (P = 


2.4) 


+0.01%,- 


-0.12% 


+0.08%,- 


-0.03% 


+0.09%, 


-0.15% 


+0.044%, -0.10% 


+0.19%,- 


-0.31% 


0.18% 


Za{P = 


2.2) 


+3.2%,- 


1.5% 


+18.1%,- 


-2.1% 


+21.3%,- 


-3.6% 


+0.17%, -0.25% 


+44.5%,- 


-8.9% 


17% 


Za{(3 = 


2.4) 


+0.6%,- 


8.0% 


+6.1%,- 


-7.3% 


+6.7%,- 


15.3% 


+4.0%, -9.2% 


+ 14.6%,- 


-33.4% 


5.9% 
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TABLE IV. Z-factors at simulation points for meson decay constants. 



p 


Z 


■TP 

V 


/7SF,L=0.8fm 


Z 


■SF,L=oc 


7 NPC 


Z 


■TP 
A 


/7SF,i=0.8fm 


r 7 SF,L = oa 

£ A 


2.187 





.81657 


0.7400(08)+"* 


0. 


.6488(24)1** 


0.4536(044) 


0. 


.83204 


yj. i i -i — i jo j 


+ 27 
-54 


0.7180(146) 


2.214 





.81923 


0.7462(07)1^ 





.6651(20)1^ 


0.4784(038) 





.83449 


7762 (341 


+24 
-49 


0.7276(127) 


2.247 





.82237 


0.7532(06)1q8 





6825(16)1?9 


0.5056(030) 





.83737 


7819(301 


+21 
-43 


0.7383(108) 


2.281 





.82548 


0.7601(06)1q7 





.6981(13)1^ 


0.5293(040) 





.84022 


0.7875(27) 


+ 19 
-38 


0.7482(092) 


2.334 





.83009 


0.7698(04)l°g 





.7187(10)li3 


0.5421(030) 





.84446 


0.7956(22) 


+ 16 
-31 


0.7619(073) 


2.416 





.83673 


0.7832(03)1^ 


0. 


,7441(07)1J 3 9 


0.5925(177) 


0. 


.85057 


0.8068(17) 


+ 12 
-23 


0.7797(052) 


2.456 





.83978 


0.7891(03)1^ 


0. 


.7544(06)1 ° 


0.6090(043) 





.85336 


0.8118(15) 


+ 10 
-20 


0.7873(045) 


2.487 





.84205 


0.7934(03)1^ 


0. 


.7617(06)+°? 


0.6420(038) 


0. 


.85546 


0.8155(14) 


+09 
-18 


0.7927(040) 


2.528 





.84496 


0.7988(02)1^ 


0. 


.7705(05)1'^ 


0.6460(038) 





.85813 


0.8201(12) 


+08 
-15 


0.7994(034) 


2.575 





.84816 


0.8046(02)1°* 


0. 


.7796(05)+* 


0.6604(034) 


0. 


.86107 


0.8250(11) 


+07 
-13 


0.8064(029) 



+ 125 

-236 
+ 109 
-203 
+093 
-171 
+080 
-143 
+063 
-111 
+045 
-076 
+038 
-064 
+034 
-056 
+029 
-047 
+024 
-039 



TABLE V. f p and in GeV with various choices of ^-factors. 



f I ^SF,Z,=0.8fmx 



0.2593(43)1!' )4 
0.2515(37)1^ 
0.2478(36)1°J 
0.2490(37)1^ 
0.2412(29)1^ 
0.2313(51)1^ 
0.2191(42)1°? 
0.2324(40)1°? 
0.2168(43)1°,? 
0.2293(36)1°,? 



f P (Z 



SF,L=oo\ 



0.2273(43)l 0! j 
0.2242(38)1^ 
0.2246(36)l O g 
0.2287(36)1^ 
0.2252(29)1^ 
0.2197(50)1^ 
0.2095(41)1°^ 
0.2232(40)1^ 
0.2091(42)1^ 
0.2222(35)1°,? 



MzY) U(z 



5F,L=0.8fm-j 



U(Z[ 



SF,L=oo^ 



0.1504(46)1^ 0.1401(65)146 

0.1446(41)1°* 0.1356(56)1^ 

0.1412(44)1°* 0.1333(56)1^ 

0.1334(36)1^ 0.1267(45)124 

0.1377(42)1°^ 0.1319(50)1^ 



2.187 
2.214 
2.247 
2.281 
2.334 
2.416 
2.456 
2.487 
2.528 
2.575 



1.017(10) 
0.966(10) 
0.917(09) 
0.896(10) 
0.829(08) 
0.734(09) 
0.674(06) 
0.652(07) 
0.612(06) 
0.574(06) 



0.2861(44) 
0.2761(38) 
0.2706(37) 
0.2704(38) 
0.2601(30) 
0.2471(54) 
0.2332(44) 
0.2467(42) 
0.2293(45) 
0.2417(37) 



0.1578(37) 
0.1601(32) 
0.1656(31) 
0.1723(33) 
0.1689(27) 
0.1744(78) 
0.1683(44) 
0.1874(40) 
0.1750(40) 
0.1872(34) 



0.1623(42) 
0.1555(37) 
0.1512(41) 
0.1423(33) 
0.1462(41) 
0.1368(39) 
0.1444(39) 
0.1358(36) 
0.1405(47) 
0.1445(53) 



0.1298(40)1° 
0.1374(40)1°. 
0.1295(37)1°: 
0.1343(47)1° 
0.1384(53)1° 



0.1254(44)1?^ 
0.1332(44)1?? 
0.1258(40)1^ 
0.1309(49)1°| 
0.1353(55)1q6 
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0.04 
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0.1340 



0.1350 
K 



0.1360 



FIG. 1. m q versus n at (3 = 2.6 on an 8 x 16 lattice. Filled symbols show measured points. The open symbol represents k c 



X 




-0.008 



FIG. 2. Effective quark mass m q (x4) at 
fitted value and error. 



= 2.6 on an 8 3 x 16 lattice. Horizontal solid lines represent fitting range, and 
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N 



0.856 



0.855 



0.854 



0.853 



0.852 



0.851 



4 8 12 16 20 24 28 32 

X, 



FIG. 3. Zv as a function of time slice at f3 = 3.125 on a 16 3 x 32 lattice. Horizontal dashed lines represent Zv at 
x 4 = T/2 = 16. 



0.870 



0.868 



N < 0.866 



oo 



0.864 



0.862 

4 8 12 16 20 24 28 32 

X 4 

FIG. 4. Values of Za as a function of time slice at /3 = 3.125 on a 16 3 x 32 lattice. Horizontal dashed lines represent Za at 
x 4 = (5/8) - T = 20. 
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FIG. 5. Time history of /a and f P at (3 — 2.4 on an 8 3 x 16 lattice. 
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FIG. 6. Time history of f A and f P at (3 = 2.6 on an 8 3 x 16 lattice. 
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FIG. 7. Histogram of fp at /3 = 2.4 on an 8 3 x 16 lattice. Intervals in :r-axis are 50. 
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0.1 



▼ Z 



SF, L=o= 



V 



0.2 

a/L 



0.3 



0.4 



FIG. 8. Volume dependence of Zv- Filled circles and solid lines are measured data and linear fits to them. Values of (3 are, 
from top to bottom, 8.0, 6.0, 4.0, 3.125, 2.8, 2.6, 2.4 and 2.2, respectively. Stars ^ F - L -°- 8fm are connected by dashed lines to 
guide eyes. Filled triangles at a/L — are Zy F ' L ~ x ' . 
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FIG. 9. Volume dependence of Za- Symbols are the same as in Fig. 8. 
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Pade fits 

SF, L=0.8 fm 



▼ Z 



SF, L=oo 
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NPC 



-v 



1 



4 



9 

FIG. 10. Comparison of Zv from tadpole improved perturbation theory and those from SF method normalized at L = 0.8 fm 
and L = oo. are determined with the conserved current. Pade fits and error are also given for Zv from the SF method. 

Vertical dashed lines represent the range we have data for decay constants. 
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FIG. 11. Comparison of Za determined by various methods. Symbols are the same as in Fig. 10. 
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FIG. 12. m q vs. /p, C ut at k c on an 8 3 x 16 lattice at (3 = 2.4. 
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FIG. 13. Z v versus m q at /3 = 2.4 on an 8 3 x 16 lattice. 
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FIG. 14. Za versus m, at (3 = 2.4 on an 8 3 x 16 lattice. 
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1, cut 



FIG. 15. Ensemble averages of /i, fv and Zy versus /i, cu t at (3 = 2.4 on an 8 x 16 lattice. 
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FIG. 16. Ensemble averages of fi, /aa and Za versus /i, cu t at (3 = 2.4 on an 8 3 x 16 lattice. 
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FIG. 17. f p versus a with various choices of Z v . Values for our improved action with Zy P and Zy V are taken from Ref. [1]. 



Values for the standard action [19] are determined with Z-factors from the conserved current. Diamond at a 
the experimental value. 
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FIG. 18. f-n versus a with various choices of Za- Symbols are the same as in Fig. 17, though f n from the standard action is 
determined with Z T A F . 
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